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■ We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ O.SSZ©) nearby dwarf galaxy, NGC 1705, in 

six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant 

circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub- 

mm by LAB OCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR 

fG measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a signicant cold dust component 

Mh within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 

^ 350 yL/m, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an 

^ additional cold dust component. 

C/^ Key words. Galaxies: ISM - Galaxies: dwarf - Galaxies: evolution - Infrared: ISM 

^"^ ■ 1 ■ Introduction As part of this effort to understand the nature of the evolu- 

^ tion of dwarf galaxies, we present HerscheQ (Pilbratt 20 1 0) pho- 

Understanding the origin and evolution of dwarf galaxies is im- tometric observation of the nearby (5.1 ± 0.6 Mpc (Tosi et al. 

portant in modern observational cosmology. Current models for 2001)) dwarf starburst galaxy, NGC 1705. The galaxy is dom- 

^ ■ galaxy formation fall into two distinct categories: construction inated optically by a massive central super star cluster (SSC) 

^ : from the "top-down" or from the "bottom-up". Which scenario NGC 1705-1 (Meurer et al. 1995), whilst studies in the mid and 

lO ; is correct? Unfortunately, the answer is still somewhat ambigu- far-IR (Cannon et al. 2006; Galametz et al. 2009) reveal the pres- 

O ; ous and left open to interpretation as both models face critical ence of two bright infrared regions flanking the central SSC, 

O . problems, in particular due to a lack of detailed, high resolution ofl'set by -250 pc from the SSC, with these ofl'-nuclear regions 

V. " studies at redshifts corresponding to the peak era of galaxy for- dominating the global IR emission of NGC 1705. This galaxy 

> mation. Bottom-up models over-predict the number and the mass provides an ideal environment for exploring the efl'ects of ongo- 

spectrum of satellites seen around galaxies such as our own and ing, massive star formation on the environment within a dwarf 

there are inconsistencies with the timescale of the build up of galaxy (Cannon et al. 2006), given its sub-solar nebular metal- 

J^ larger galaxies (e.g., Prantzos & Silk 1998). A major source of Hdty (Z -35% Z©; (Lee et al. 2004) and large reservoir of gas 

this ambiguity may arise from the fact that accurate measure- (Meurer et al. 1998), as we can trace the efl'ects of the SSC on 

ments of the interstellar medium (ISM) properties and star for- the surrounding interstellar medium, and in particular, can char- 

mation histories of dwarf galaxies have been lacking for the most acterize the nature of dust and PAH emission, 
part. To address this crucial issue, we have begun to focus our 
study of galaxy evolution onto analogues of these earliest galaxy 

building blocks right on our door step - nearby, low metallicity 2. Observations 
dwarf galaxies - where we can begin to discern the evolutionary 

processes at work in these objects at very high spatial resolution. NGC 1705 was observed as part of the Dwarf Galaxy Survey 
Such objects are the most common type of galaxy in the cur- programme (PI. S. Madden), a Guaranteed Time (GT) key pro- 
rent epoch, making up 95% of the Local Group galaxies alone gramwiththeobjectiveof mapping the dust and gas in 51 nearby 

(Mateo 1998). 

^ Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 



X 



email: b.ohalloran@imperial.ac.uk tant participation from NASA. 



O'Halloran et al.: Herschel observations of NGC 1705 



-53" 22' 30" - 



-53" 22' 30" - 



3.6 /2m 


8 /xm 




24 /i.m 


>- 


ssc 

1 






% 


1 

D1 




• •^ 


70 fiui 


100 fim 


_ 


160 /i.m 


*- 


•* 


- 


* 










" ® 




250 fim 


350 fim 


■ ■ ^ 




■ ■ ' 


00 /xm 


^^ 








^^^1 






^^^2 










" 








1 


J 



04"^ 54™ 20^ 15^ 10^ 04^ 54™ 20"= 15^ 10^ 05^ °4 54 20 15 

"J2000 "J2000 



Fig. 1. IRAC/MIPS (top) PACS (middle) and SPIRE (bottom) maps of NGC 1705. The top three panels are the Spitzer IRAC and 
MIPS maps (1-r) at 3.6, 8 and 24 yum, the middle panels are the PACS maps (1-r) at 70, 100 and 160 yum, while the bottom three 
panels are the maps at 250, 350 and 500 yum respectively. The central SSC and the eastern (Dl) and western IR components are 
noted on the IRAC 8 yum map, whilst the beam FWHM values for the PACS and SPIRE maps are also illustrated. 



dwarf galaxies, sampling a broad metallicity range of 1/50 to 1/3 

Zo. 



2A. SPIRE observations and data processing 

The galaxy was observed by SPIRE (Griffin et al. 2010) at 250, 
350 and 500 yum for a total of 733 seconds, in scan-map mode 
with scanning rate 30"/sec, with the final map covering roughly 
16x16 arcmin. The measured 1 cr noise level are 5, 6 and 7 mJy 
beam"^ at 250, 350 and 500 yum respectively; the noise levels in 
the images are dominated by confusion. The data were processed 
using the HIPE pipeline (see Pohlen et al. (2010) for a detailed 
description, Swinyard et al. (2010) for calibration accuracy and 
Bendo et al. (2010b) for details on the destriper). The pipeline 
produces maps with a pixel size of 6.0, 10.0 and 14.0" at 250, 
350 and 500 yum respectively. The ICC has released some interim 
small correction factors to improve the preliminary calibration. 
All flux values derived using the current standard calibration file 
for the flux conversion, are multiplied by 1 .02, 1 .05, and 0.94, for 
the 250yum, 350yum, and 500yum maps, giving final global fluxes 
of 0.85 ± 0.13, 0.38 ± 0.06 and 0.26 ± 0.04, respectively, using 
an aperture of 72" as per Galametz et al. (2009). The uncertainty 
in the flux calibration is of the order of 15%. 



2.2. PACS data processing 

The galaxy was observed by PACS (Poglitcsh et al. 2010) at 70, 
100 and 160 yum in scan map mode for a total of 1504 seconds, 
with the final map covering roughly 7x6 arcmin. The PACS 
data were reduced starting from the Level product using the 



HIPE version 3.0 data reduction software and the HIPE respon- 
sivity calibration file version 1, corrected by factors provided 
by the ICC. We applied the basic steps of the standard reduc- 
tion pipeline to perform the data reduction to the Level 1 , where 
we mask the bad and saturated pixels, apply a flat field correc- 
tion and perform a multi-resolution median transform (MMT) 
deglitching and apply a second order deglitching procedure to 
the data. We mask the data contributing to the bright structures 
in the data cube and perform polynomial fits on half scan legs to 
subtract the baselines. This process removes most of the drifts 
of the maps. Two dimensional maps are finally constructed us- 
ing the MadMap procedure of HIPE. Striping is visible in a sin- 
gle scan but is mostly removed by the cross scan. The pipeline 
produces maps with a pixel size of 3.2, 3.2 and 6.4" at 70, 100 
and 160 yum respectively. In this version of the responsivity cal- 
ibration file, an error still exists in the flux scale. The final re- 
duced global fluxes have to be scaled down by 1.05, 1.09 and 
1.29 for the blue (PACS 70), green (PACS 100) and red (PACS 
160) bands, giving 1.05 ± 0.12, 1.22 ± 0.13 and 1.18 ± 0.12 
Jy, respectively, using an aperture of 72" as per Galametz et al. 
(2009). The uncertainty in the flux calibration is of the order of 
10%. 



3. Results & Discussion 

3.1. IVIorphology of NGC 1705 

The reduced SPIRE and PACS maps of NGC 1705 are pre- 
sented in Figure 1. To perform aperture photometry on indi- 
vidual sources with NGC 1705 for both the SPIRE and PACS 
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Fig. 2. Global spectral energy distribution for NGC 1705, with 
a final cold dust temperature of 5.8 K andyS =1. The grey curve 
denotes the combined dust model, the green denotes the stellar 
component, the red denotes the warm dust/PAH components and 
the blue denotes the cold dust component. The symbols refer to 
the Herschel and Spitzer fluxes used to fit the models, in con- 
junction with the 850 yum LABOCA flux. 



images, we used the IRAF package apphot, in which the back- 
ground is estimated using annuli just outside the boundaries of 
our galaxy. To determine the flux densities at MIPS bands, as a 
comparison to the PACS fluxes, the observations are convolved 
and re-gridded to a common resolution (FWHM MIPS 160 yum: 
40"). We use convolution kernels (Gordon et al. 2008; Bendo et 
al. 2010) which convert a higher resolution PACS point- spread 
function (PSF) to lower resolution MIPS PSFs using Fourier 
transforms. 

As can be seen in both sets of images, the morphology of 
NGC 1705 is quite complex - the far-IR emission is domi- 
nated by the ofl'-center sources previously noted by Cannon et 
al. (2006) and Galametz et al. (2009). As might be expected, 
the two IR regions are clearly resolved at the shortest PACS 
wavelengths, with no emission detected from the central SSC 
which indicates very little in the way of significant dust emission 
within the central SSC itself. As we move into the SPIRE bands 
(beyond 250 yum), the individual components are no longer re- 
solved. As noted in Galametz et al. (2009), the off'-center bright 
IR fluxes do not have stellar counterparts. There is no indication 
of the presence of the ofl'set peak to the west of the galaxy in the 
PACS or SPIRE that had been previously noted by Galametz et 
al. (2009), nor are there counterparts to the faint 24 yum sources 
also noted. Of note is the fact that the easternmost IR source 
(Dl, in the nomenclature of Cannon et al. (2006)) is the bright- 
est across the difl'erent resolved wavelength maps, and is of the 
order of 1.15-1.45 times brighter than the western source, with 
the flux ratio peaking at 100 yum - previous IRS spectroscopy 
by Cannon et al. (2006) confirms the strongest PAH emission in 
this region, which in conjunction with the Herschel data is sug- 
gestive of a greater concentration/contribution to the dust mass 
within this region versus both the westernmost region and the 
SSC. 



3.2. Spectral energy distribution of NGC 1705 

3.2.1 . Cold dust component scenario 

To quantify the contributions to the IR luminosity from various 
dust populations within NGC 1705, we attempted to fit modified 
blackbody functions to the SPIRE and PACS measurements in 
an eff'ort to obtain dust component luminosities, masses and dust 
temperatures. On fitting, it became quickly apparent that fitting 
a single blackbody to the SPIRE and PACS fluxes was not possi- 
ble, as the SED diverged from a single temperature fit longward 
of approximately 350 yum with the 350 and 500 yum fluxes higher 
than the fit, a trend that was confirmed with the addition of the 
LABOCA 850 yum flux from Galametz et al. (2009). Even the 
addition of a second modified blackbody function to the fit was 
unable to fit to a reasonable degree the 350 and 500 yum fluxes in 
conjunction with the LABOCA flux - we tried to fit a number of 
modified blackbodies with a range of realistic emissivity values 
between 1 and 2, and were unable to obtain a satisfactory fit that 
suitably fit the three longest wavelength fluxes. 

In light of this, we instead constructed a global spectral en- 
ergy distribution for NGC 1705, which in addition to the SPIRE 
and PACS fluxes includes global fluxes from LABOCA and 
the Spitzer IRAC and MIPS fluxes from Cannon et al. (2006) 
to better constrain the mid-IR contribution to the SED, and 
this global SED is presented in Figure 2. We used a realistic 
dust radiation model as far as dust grain properties were con- 
cerned, a simpler version of the Galliano et al. (2008) model. The 
sources of emission are the dust grains (molecules like PAHs, 
silicates, graphites) and old stars, while the dust composition 
and size distribution adopted were those of Zubko et al (2004), 
with abundances assumed to be solar. The dust optical proper- 
ties were taken from Draine and Li (2007) for PAHs, Laor & 
Draine (1993) for graphites and Weingartner & Draine (2001) 
for silicates. The interstellar radiation field spectral shape was 
assumed to be that of the Galactic diff'use ISM (Mathis et al 
1983). The Dale et al. (2001) prescription was used to link the 
dust mass to the integrated radiation density it was exposed to 
- see Galametz et al (2009) for a more complete description of 
the free parameters of the modelling, among those the total dust 
mass (Mdust) and the PAH-to-dust mass (fpAH) nomalized to the 
Galactic value. To account for the observed excess at submm 
wavelengths, an additional, independent, modified blackbody 
was added to the model. 

We then proceeded on two fronts in terms of fitting the model 
to the fluxes - a) letting /3 vary and keeping the temperature 
fixed and b) keeping JS fixed, and letting the temperature vary. 
However, we found that whcnjS was allowed to vary with a tem- 
perature of the cold dust fixed (lOK), the SED fit had a tendancy 
to find beta values less than 0.5. This unrealistic value ofjS (nor- 
mally with values between 1 and 2) (Bolatto et al., 2000) prob- 
ably reflects instead, a distribution of dust temperatures. As a 
result of this, we decided to fix the value of /? to 1 and let the 
temperature of the cold dust component vary instead. 

The addition of such a dust component (with a temperature 
of < 10 K) within the global SED of NGC 1705 has been ten- 
tatively suggested previously by Galametz et al. (2009), on the 
basis of the addition of the 850 yum LABOCA point, with this ad- 
ditional dust component containing up to 70% of the total dust 
content of the system. There have been significant claims of prior 
detections of very cold dust components within dwarf galaxies, 
(GalHano et al. (2003, 2005) and Galametz et al. (2009)), but the 
lack of far-IR and sub-mm fluxes to constrain SEDs sufficiently 
well has left this open to debate. If present, very cold dust can 
only be found deep in the interior of clouds; it must be protected 



O'Halloran et al.: Herschel observations of NGC 1705 



from stellar light and shielded by an optical depth well above 
10 mag. The PACS and SPIRE images provide us with an addi- 
tional constraint regarding the cold dust, namely it must be lo- 
cated where we see the emission, e.g. the eastern condensation, 
which strongly is suggestive of high extinction in those loca- 
tions. With the addition of the three SPIRE photometry points, 
in particular the 350 and 500 fim points, the claim for an addi- 
tional dust component is strengthened. 

Our best-fit modeling is consistent with the need for an addi- 
tional cold dust component, and leads to a cold dust mass of 4.0 
X 10^ M© with a temperature of the cold dust of 5.8 K (withyS = 
1). For the warmer dust component, we estimate a dust mass of 
2.1 X 10^ Mo, giving a total dust mass of 4.2 x 10^ M©. The plot- 
ted fit in Figure 2 is indeed the best available fit - we attempted 
to fit a range of dust temperatures to satisfy a reasonable fit to the 
submillimeter excess. However, it was not until the dust temper- 
atures reached less than lOK that we began to obtain acceptable 
fits - we considered a cold dust component fit with a temperature 
of lOK and a reduced chi squared value of 1.1 to be the lower 
limit in terms of an acceptable fit. In terms of obtaining a reason- 
able lower limit for the dust mass, fitting a cold dust component 
of 10 K leads to an acceptable reduced chi squared value of 1.1, 
giving us a cold dust mass of 1.1 x 10^ M© and a total dust mass 
of l.SxlO^M©. 

The multiband Herschel observations have allowed us to 
constrain significantly the SED compared to the previous work 
of Galametz et al. (2009), where the only sub-mm constraint was 
the LAB OCA flux, allowing us to refine the dust component 
masses and by extension, obtain new values for the PAH/dust 
and dust to gas ratios. For the MpAH/^dust, we get ~ 2.3 x 10""^ - 
much lower than the Galactic value, and consistent with the PAH 
emission deficit at very low metallicity (Madden 2002, Galliano 
et al. 2003, O'Halloran et al. 2006, 2008). Correspondingly for 
the dust to gas ratio, we get ~ 1.0 x 10"^, using the HI mass 
from Galametz et al. (2009). The derived dust/gas ratio is rather 
lower than expected for a low metallicity dwarf galaxy - we 
could simply be predicting too much dust within our SED mod- 
eling, as a result of the use of graphite grains. Serra Diaz-Cano 
& Jones (2008) warn against the use of graphites in construction 
of SED modeling, suggesting instead that amorphous carbons 
could be more relevant to represent the interstellar carbon dust, 
since they possess a higher emissivity at submm wavelengths 
and thus require less mass to account for the same emission. 
Substituting in amorphous carbons to replace the graphite grains 
in our model (using values from Rouleau & Martin (1991)) and 
using the same parameters of yS and T as before, we obtain re- 
duced warmer (2.86 x 10"^ M©) and cold dust masses (9.8 x 10"^ 
M©), with a total dust mass of 1.2 x 10^ M©. These reduced dust 
masses have the knock-on eff'ect of a higher dust-to-gas ratio (~ 
2.3 X 10"^), a value that is much closer to that expected for such 
a low metallicity galaxy. Both D/G values (from using graphites 
and amorphous carbons) are still at the lower end of the chem- 
ical evolution model used in Galliano, Dwek & Chanial (2008) 
- nonetheless, these are values which seem consistent with an 
object that while possessing a low metallicity, does contain sig- 
nificant amounts of warm and cold dust. 



3.2.2. Alternatives to the cold dust component scenario 

We must also consider the possibility that the sub-mm excess de- 
tected in NGC 1705 is not as a result of an extra dust component. 
One hypothesis arguing against such a conclusion is that the sub- 
mm excess could originate instead from hot (~100 K) dust with 



a dust emissivity index =1 and the temperature fluctuations of 
very small grains (Lisenfeld et al. 2002). 

The excess could also be explained through a change in the 
emissivity of the cold dust grains. For an example of how the 
emissivity may change as we move into difl'ering temperature 
regimes, Dupac et al. (2003) suggest that JS decreases with in- 
creasing temperatures, from ~2 in cold (T ~ 1 1-20 K) regions to 
0.8 to 1.6 in warmer regions (T ~ 35 - 80 K). In contrast, Paradis 
et al. (2009) showed that the spectral shape of emissivity are al- 
ways steeper in the FIR (A< 600 yum) and flatten in the submm 
and mm regions. In regions where dust is significantly colder in 
the molecular phase than in the atomic phase, an increase in the 
emissivity by a factor of ^3 was detected only in the FIR, whilst 
the emissivity for the dust in the atomic and molecular phases 
become comparable again in the submm and mm wavelength 
range. The observed break in the emissivity spectrum is in qual- 
itative agreement with the dust emission model of Meny et al. 
(2007), which invokes quantum eff'ects in amorphous solids to 
explain the flatness of the observed submm emission spectrum 
and also produces a break in the emissivity slope around 600 
yum. 

However, one must exercise caution in adopting such an in- 
terpretation. Flux uncertainties, especially in the Ray leigh- Jeans 
regime, can afl'ect the results for the SED fits as far as tempera- 
ture and emissivity are concerned, as fitting fluxes near the SED 
peak produces inaccurate temperature and dust spectral index es- 
timates. In addition, line-of- sight temperature (and density) vari- 
ations can also afl'ect the SED fitting (Shetty et al. (2009a,b)). 
Longer wavelength fluxes in the Ray leigh- Jeans part of the spec- 
trum (> 600 yum) may more accurately recover the spectral in- 
dex, but both methods are very sensitive to noise (Shetty et al. 
(2009a,b)). 

An additional alternative to a cold dust component is spin- 
ning dust where the rotational dynamics of very small interstel- 
lar grains can explain the 10 - 100 GHz component of the difl'use 
Galactic background via electric dipole rotational emission un- 
der normal interstellar conditions (Draine & Lazarian 1998a,b). 
However, Jones (2009) notes that observations by Dickinson et 
al. (2006), which searched for a microwave emission excess 
in an H II region, puts an upper limit on the dust emission at 
31 GHz and appears to be inconsistent with the spinning dust 
model. Given the unsatisfactory evidence arguing against the use 
of amorphous carbon to explain the sub-mm excess, we conclude 
that the most likely scenario to explain the observed excess is 
that of an additional cold dust component. 
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